Study of the 2d Ising Model with Mixed Perturbation 

p. Grinza"*and A. Rago"^ 

Dipartimento di Fisica Teorica deU'Universita di Torino and 
Istituto Nazionale di Fisica Nucleare, Sezione di Torino 
via P.Giuria 1, 1-10125 Torino, Italy 

February 1, 2008 



Abstract 

We study the thermodynamical observables of the 2d Ising model in the neigh- 
borhood of the magnetic axis by means of numerical diagonalization of the transfer 
matrix. In particular, we estimate the leading order corrections to the Zamolod- 
chikov mass spectrum and find evidence of non-vanishing contributions due to the 
stress-energy tensor. 
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Introduction 



In the past the 2d Ising model has been the subject of both analytical and numerical 
study. After the original solution due to Onsager much work has been done to improve 
our knowledge of the model. 

It has been shown that at criticality the model can be described by a Minimal Unitary 
Conformal Field Theory 0. Both magnetic and thermal perturbations of this model 
have been investigated and it has been found that these are the only two integrable 
perturbations of the model. The integrability of the model with a magnetic perturbation 
has been exploited by Zamolodchikov in order to obtain the exact S-matrix and mass 
spectrum of the theory 0. Further, the exact knowledge of the S*- matrix has been utilized 
to calculate the first few terms in the spectral expansion of the two point correlation 
functions via the exact calculation of the form factors 

On the other hand, an infrared safe short distance expansion (IRS) for the correlators 
has been proposed in (see |^ for a comparison of the expansions to Monte Carlo 

simulations). An interesting feature of the latter approach is that, in principle, it does 
not require the integrability of the model. 

Apart from the analytic results there exists an extensive numerical investigation of 
the model. In particular, the transfer matrix technique has been employed in a large 
set of investigations: the study of the magnetic perturbation of the model P, even at 
finite temperature 0], the critical equation of state [|l^ and a classification of irrelevant 
operators which enters in the observables of the theory ||TT],|r^. 

However, the knowledge of the mixed perturbation of the model is still limited. The 
most important contributions to the study of this regime are due to McCoy and Wu |T3], 



Delfino et al. 14 and to Zamolodchikov and Fonseca [15 



This work explores this latter regime, in fact, our purpose is to investigate the mixed 
perturbation model in the neighborhood of the magnetic axis. We apply the method 
developed by [§ to this case. There are two key points of this method: the high precision 
numerical determination of the eigenvalues of the transfer matrix as a function of the 
couplings and the analysis of the data by means of scaling functions obtained from a CFT 
approach. 

The original contributions of this work can be summarized in three points. Firstly, 
we verify that the CFT inspired scaling functions are in perfect agreement with our data, 
moreover we are able to estimate the leading order corrections induced by the presence 
of the thermal perturbation. 

Secondly, we have been able to predict the existence of a term due to the stress- 
energy tensor in the free energy and determine its amplitude. This is an interesting result 
because the contributions of the stress-energy tensor have been claimed to be zero if only 
integrable perturbations are involved. However, this is issue is still a matter of study. 

Finally we give an estimate of the corrections to the Zamolodchikov mass spectrum 
and we find perfect agreement (up to some acceptable error) with Delfino et al. [|14|. 

This work is organized as follows. In Section |I] we review the standard definitions used 
in the Ising model and we define the normalizations and conventions we use. Section 
is devoted to the explanation of the transfer matrix technique. The computation of the 
scaling functions is illustrated in Section ^. Our results are presented in Section ||, and 
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finally, we draw our conclusions in Section ^. In Appendix |^ we summarize all of the 
known results on the amplitudes of the scaling functions, in Appendix |B] we present an 
example of the obtained scaling functions. 



1 The Ising Model 

In this section we review the existing results of the 2d Ising model. In particular, we shall 
present the model both in its lattice form and in its continuum formulation as a field 
theory. We also define the observables that we shall use in the following work. 



1.1 The Lattice Model 

The Ising model in a magnetic field at an arbitrary temperature is defined by the partition 
function 

Z{l3,hg) = e^^<-->''"'''"+''^^"''", (1) 

(Ti=±l 

where the spin variable cx^ takes the values ±1; the notation {n,m) represents nearest 
neighbor sites on the lattice; the sites are labeled by n = (no, rii) and the two sizes of the 
square lattice are Lq and Li (they are taken to be different because our transfer matrix 
calculations will treat the two directions asymmetrically); the total number of sites of the 
lattice will be denoted as A^ = LqLi. 

The coupling (3 is the inverse of the temperature, while the magnetic perturbation is 
introduced by the coupling h£ = HjS, where H is the magnetic coupling. This model 
undergoes a second order phase transition when he = and P reaches its critical value jSc 

Pc = ^ log( V2 + 1) = 0.4406868 .... (2) 

Now we define the observables which we shall consider in the following work. 
• The Magnetization per site is defined as 



The Magnetic Susceptibility is defined as 



X(AM = ^^- (4) 



The Free Energy is defined as 



/(/3,M = llogZ(/3,M. (5) 
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The free energy is composed of a "bulk" term fb{(3, hi), which is an analytic function 
of hi and /3, and a "singular" term encoding the relevant information about the 
theory in the neighborhood of the critical point. Prom the exact solution of the 
lattice model &t — and /3 = /3c we can compute the value of /^(O, 0) 

2(7 -[ 

/b = — + - log 2 = 0.9296953982 . . . , (6) 

TT 2 

where G is the Catalan constant. 

• The Internal Energy density is defined as 

mhe) = ^{J2^nam). (7) 

(n,m) 

As for the free energy, we have both a bulk and a singular part. If we define £'b(/3, /i^) 
as the bulk contribution, we can obtain the value of = Eb{0, 0) = (by using 
Kramers- Wannier duality or CFT techniques) and we can also define (for future 
convenience) 

E{f3,he)=E{P,he)-eb. (8) 

• Time slice correlation functions. 

Wc can define the zero momentum projections of the two-point correlation func- 
tions {a{r)a{0)) and (e(r)e(O)) (they are also called time slice correlators). The 
magnetization of a row of the lattice (time slice) is given by 

^no ^ J^^Cr{no,ni) (9) 

hence, the correlation function between time slices is given by 

Glirir) = ^ [{Sn,S^,+n.) - (5„,)^] (10) 
no 

where label indicates that it is the zero momentum projection of the original 
correlator. 

1.2 The Ising Field Theory 

The 2d Ising model near the phase transition can be described via a minimal unitary 
Conformal Field Theory (with central charge c = 1/2) perturbed with the energy and 
magnetization densities e{x) and a{x) 

A = AcFT + h J d'^xa{x)+T J d'^x e{x). (11) 

Where Ac ft is the action of the model at the critical temperature without external mag- 
netic field. 
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The coupling constant h represents the magnetic perturbation of the model and, as stated 
before, it is the continuum version of the previously defined coupling he- The other cou- 
pling constant r represents the thermal perturbation and, near criticality, it is proportional 
to the reduced temperature t 

'=\^ (12) 
which we use in the rest of the paper. 

In the following Section we give a brief report of the known results on the field theory 
description of the model. 

1.2.1 The Critical Theory {t = 0, /i^ = 0) 

The Ising model is the lowest model of the so-called "Minimal Unitary" series of conformal 
theories whose central charge is given by 

c, = l--^, p=3,.... (13) 

The pecuharity of these models is that they possess a finite set of primary fields; as a 
consequence the whole space of local operators of the theory can be built by applying the 
generators of the Virasoro algebra to the primary fields. The operators obtained in this 
way are called secondary fields or descendants. Following this route, one is led to organize 
the operator content of the theory in conformal famihes, i.e. the sets of descendants of 
each primary field. 

The operator spectrum of the Ising model consists of three primary operators 

• Identity ^ A/ = 

• Magnetization =^ A^. = 1/16 

• Energy =^ A^ = 1/2 

where Xi — 2Aj is the scaling dimension of each operator. 

We assume for the field a{x) and e(x) the usual CFT normalization 

(a(x)a(O)) = \x\-^/\ {e{x)e{0)) = \x\-^ (14) 

so they have scaling dimensions 1/16 and 1/2 respectively; the parameters h and t have 
dimensions 15/8 and 1. 

Hence there are three conformal famihes that descend from these operators and the general 
expression for the descendants 0[<^] is 

= L_k^...L_k^L_^^...L_m^(f) (15) 

where 

s p 

^ki = n; ^mi = n, (16) 

i=l i=l 
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the operator (f> is one of the primary fields of the theory. 

The scahng dimension and the conformal spin of the operators are given by 

Xo = Ao + Ao =n + n + A^ + A^ ^^^^ 

So = n — n. 

Among secondary fields, the quasi-primary fields play a special role. A descendant field 
Q is called quasi-primary when 

• L,\Q) = 

• IQ) is not a null vector. 

It will be shown in Section |^ that quasi-primary fields can be used as building blocks for 
the construction of an effective Hamiltonian for the lattice model near criticality. 

1.2.2 The Ising Field Theory with Magnetic Perturbation {t = 0,he ^ 0) 

If (3 is fixed to its critical value Pc = |log(l + V2) and the magnetic field is switched 
on, the field theory is still integrable, i.e. it possesses an infinite number of integrals of 
motion. This implies the exact knowledge of the S'-matrix and the mass spectrum. The 



latter result, due to Zamolodchikov [[l^], suggests that the model could be described by a 



scattering theory with a spectrum of eight self-conjugated particles with masses 



7712 


= 2mi 


cos 




1.6180339887..) 7771, 


ma 


= 27771 


cos 


TT 

30 ~ 


(1.9890437907..) 7771, 


7774 


= 27772 


cos 


7n 

30 ~ 


(2.4048671724..) 7771, 


77i5 


= 27712 


COS 


277 

15 ~ 


(2.9562952015..) mi. 


7776 


= 2r7l2 


COS 


TT 

30 ~ 


(3.2183404585..) mi. 


7717 


= 4m2 


COS 


TT 

— COS 

5 


^ = (3.8911568233..) mi 


7778 


= Am2 


COS 


TT 

— COS 

5 


^ = (4.7833861168..) mi 



:i8) 



where mi = mi{h) is the fundamental mass of the theory, and is given by 

mi{h)=ChT5. (19) 



The numerical value of C was computed by Fateev [17 



4sin^r(i) /47r2r(|)r2(l|)- ^ ^ 

C = 5 V5^ n\ ,^ 3 \ = 4.40490858 .... (20) 



r(i)r(^) V nm,. 

The vacuum expectation values of energy and magnetization can be parametrized as 

(e) = A, /7«/l^ {a) = A„ h'^'^ (21) 

where the amplitudes and A^- can be computed exactly |jl8 



A, = 2.00314..., A^ = 1.27758227.... (22) 
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1.2.3 Non-integrable Perturbation of Ising Model (t 0,he 0) 

The non-integrable perturbation of the 2d Ising modelQ was treated in [|l^ in the frame- 
work of the Form Factors Perturbation theory. 

The starting point of FFPT is to consider the mixed perturbation as the perturbation of 
an integrable QFT. 

The action (Eq. [Tl|) describes a one parameter family of theories labeled by the adi- 
mensional scaling variable ^ 

e = -4- (23) 

/lis 

In this framework it is possible to calculate (at first order) the corrections to the mass 
spectrum and vacuum energy of the integrable theory. Two limits have been investigated 
14| : the magnetic perturbation of a free massive Majorana fermion oo) and the 



m 



thermal perturbation of the integrable magnetic perturbation of the critical Ising model 

(e-0). 



The detailed discussion of shows that in the former case oo), a straightforward 
application of the method is not an easy task. In fact, if we are in the high temperature 
regime of the model, the corrections at first order vanish. A second order calculation is 
required but is a non-trivial computation. However, in the low temperature regime, the 
divergence remains at first order, and is the signal kink's confinement. 

The latter case — > 0), however, is more tractable and gives rise to some quantitative 
predictions of the mass spectrum of the model. The results are listed below 

^^'•"'^ -0.0558... m?. 



6 mi 
6m I 
6mi 



0.8616..., (24) 
1.5082.... 



The non integrable perturbation gives rise to an important quantitative difference on the 
masses above the threshold. In fact, the integrability of the pure magnetic perturbation 
prevents the creation of new particles at energies above the threshold (and from any 
other inelastic process). The explicit breakdown of integrability implies that particles 
above threshold become unstable and are expected to decay. This new feature manifests 
itself in the corrections to the masses, namely they develop an imaginary part. This effect 
has been seen explicitly in the case of the first mass above the threshold m^. The first 
order correction is real and given by 

6mA , , 

— ^ ~ 1.1460... 25 

while, at second order, one expects a non- zero value of Im (it was not computed in |T4 
because the authors' analysis only covers the first order contribution). 



"'^See, e.g. for a detailed discussion of the model with thermal perturbation only (t 7^ 0, /i^ = 0). 
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1.3 Relation Between Lattice and Continuum Operators 



It is useful to define, for future convenience, the relations between the lattice and contin- 
uum definitions of the energy and magnetization operators. Near the critical point, the 
simplest choices for the lattice operators are 

• Spin operator 

= (26) 

where the index I indicates that it is a lattice discretization of the continuum oper- 
ator. The magnetization per site is defined as 



N 



(27) 



Energy operator 



y&{n.n.x} 



(28) 



where the sum runs over the nearest neighbor sites of x and e{, is the bulk term. 
The energy per site is defined as 



(29) 



Off-critical corrections to scaling of these operators will be discussed later. 



1.4 Converting between Lattice and Continuum Units 

In order to fix the conversions between lattice and continuum units, we shall follow the 
careful discussion of One can write the lattice versions of a and e as 

ai = f^{t,h)cr + f,{t,h)^, ^^^^ 
e« = 9oit, h)e + gi{t,h)(t), 

where Jq , fi, Qq, Qi are suitable functions of t and h (which also depend on the parity prop- 
erties of the operator), while the operators 0i are all other fields (relevant and irrelevant) 
of the theory respecting the symmetries of the lattice. We also have 

he = bo{t,h)h (31) 

which is the relation between the lattice coupling constant and the continuum magnetic 
field h; bo{t, h) is an even function of h. 

At first order in t and i.e. near criticality when t — >■ and /;,—>■ 0, we have 

ai = R^a, ei = R^e, hi = Rhh (32) 
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where the constants Ra, Re, Rh are defined as 

R^= ljmf^{t,h) 

= tfeo^°^^'^^ (33) 
Rh = hm boit, h). 

The previous normahzations were fixed in ^ by comparison with the exphcit expression of 
the spin-spin and energy-energy critical correlation functions on the lattice. The numerical 
results are 

R„ = 0.83868 . . . 

i?^ = i = 0.31831 . . . (34) 
vr 

= = 1.1923.... 



2 Transfer Matrix Technique 



We face the problem of computing mass spectrum and observables by numerical diagonal- 
ization of the transfer matrix. This technique, introduced in 1941 by Kramers and Wan- 
nier pO[ , was extensively used by Baxter to obtain analytic solutions of certain statistical 



mechanical models For further details and discussions about numerical implementa- 
tions of the transfer matrix see, e.g. . 



The basic idea is to rewrite the Boltzmann weight by means of the transfer matrix T(uj, Uj 



T{Uno,Uno+l) = ViUnoY^^ f/(«no , ^no+l) ^(^no+l) 



1/2 



(35) 



with 



U{uno,Uno+i) = exp /3 ^ 

V ni=l 
/ Li 

^(■"no) = exp I P ^ 0-(no,ni)0"(no,ni+l) + (np ,ni) 



(36) 



ni=l 



ni 



where = {o-{no,i),o'{no,2), ■ ■ ■ 5cr(no,Li)) IS the spin configuration at the row (time-slice) 
Uq. The previous position implies that the partition function becomes 

= ^ T{ui,U2) T{U2,U3) ■■■ T{ulo,Ui) 



(37) 



where T is a positive and symmetric 2^^ x 2^^ matrix, whose (real) eigenvalues are the 
A,. 
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2.1 Observables in Transfer Matrix Formalism 

The numerical computation of eigenvalues and eigenvectors of the transfer matrix enables 
us to compute all the observables we need, provided that we specify the values of hi, t 
and Li. We can derive all the observables from the partition function Z{j3, hi) (|37|). For 
the derivation see iSll. 



Free energy 



f{P,hi) = -J— log Z {(3, hi). (38) 



This expression simplifies in the limit Lq ^ oo, in fact the leading contribution is 
due to the maximum eigenvalue Xmax 

f{(3,hi) ~ ^logA^,,. (39) 

Magnetization 

, ^ tr5T^« 

(ai) = -— 40 

where S{uno,Uni) = S{uno, ^ini)o'(„o,i). In the limit Lq oo we have 

Wi) = m\o) (41) 

where |0) is the eigenvector associated to the Xmax eigenvalue. 
Energy 

The case of the internal energy is analogous similar to that of magnetization. In the 
limit Lq — > oo we have 

(q) = {0\E\0) (42) 
where the matrix E{uno,Uni) is given by 

E{Uno, Uni) = S{Uno, Uni)o-(^no,l)0'{no,2)- (43) 

Correlation functions and mass spectrum 

In the limit Lq —>■ oo the time slice correlation function is defined as 

{So St) = 5^exp(-m, \t\){0\S\i) {z\S\0) (44) 

i 

with S = S{Uno,Un,) E„i ^(no.ni)- 

The mass spectrum is given by 

mi = -log(^] (45) 



A, 



where the eigenvalues are organized in decreasing order of magnitude Xmax = Xq > 
Ai > ■ ■ ■ > Aj > . . . and \i) are the normalized eigenvectors of T. 
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3 Scaling Functions 



In order to study the collected data obtained from the TM, we need to know the behavior 
of the measured operator as a function of the perturbation variables. The fundamental 
ingredient of this construction is the knowledge of the whole spectrum of operators of the 
theory, including the OPE between them. 

Hence, the operator content of the 2d Ising model at the critical point previously discussed, 
enables us to build an effective Hamiltonian for the perturbed model. As discussed in 
detail in ||^, the aim of this effective Hamiltonian is not to describe the model at a scale 
comparable with the lattice spacing; instead it has to be considered as the Hamiltonian 
describing the model after a suitable number of Renormalization Group transformations, 
i.e. at a scale that is larger with respect to the lattice spacing. 

3.1 Lattice Construction of the Ising Model Via CFT Operators 

The main idea then is to use the whole spectrum of conformal operators, defined on the 
continuum, to describe the corrections to scaling (due to the lattice) in the observables of 
the model. In order to build this Hamiltonian explicitly, we have to take into account, in 
principle, the following ingredients: 

• Symmetries of the model 

Unlike the case of the model at the critical point, which exhibits two exact sym- 
metries (Z2 and duality), the presence of the magnetic field explicitly breaks all of 
them. Hence, in this case, there are no constraints coming from symmetries (it is 
useful to remember that this argument, in the critical case, selects only the fields 
belonging to the conformal family of the identity) . 

• Symmetries of the lattice 

It is crucial to define the geometry of lattice we are using in the transfer matrix 
calculation. In the following we consider a square lattice. This means that the 
rotational symmetry of the CFT is broken down to the dihedral subgroup D4 and 
also operators with spin are allowed. Hence, the residual symmetry group (rotations 
of integer multiples of 7r/2) implies that only operators with spin j = 4k, k & N, 
can appear on the lattice (see |^ for a detailed discussion). 

• Lattice continuum relations 

Lattice operators are defined in terms of continuum operators as follows 



where / and g are functions of the reduced temperature t, and they are even func- 
tions of the magnetization he. Furthermore, the operators appearing in the previous 
expressions have to be compatible with the described symmetries. 

Now we are able to write the following lattice (effective) Hamiltonian 



ai = f^ihe, t)a + hf^ih, t)e + /f (/i,, t)ai + /i,/;(/i,, t)^ + h,flr]i 
Q = gl{hi,, t)e + hg^ihi, t)a + heg^{h(,, t)ai + gl{hi, t)ti + g]rii 



(46) 



Hiat = He FT + hia + te + Ui^!i 



(47) 
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where are the quasi-primary fields belonging to the whole set of conformal families of 
the theory with spin j = 4k, A; G N. 

The least irrelevant fields which enter the expression of Hiat are built starting from the 
following quasi-primary fields of the family of the identity 



Ql = L_2l 



5 



(48) 



where the notation Ql is used to denote with rj the conformal family and with n the 
level of descent. The same can be done for all the other families (for a list of low-lying 



quasi-primary states, up to level 10, see |TT 
All the results are reported in Table |I]. 



Identity 
Energy 
Spin 



Spin-0 Sector 


Spin-4 Sector 


RG 
Eigenvalue 


QlQl = Tf 


gi + = T2 + f 2 


-2 


Q4Q4 


^56^2 + QlQ\ 


-6 


Q4Q4 




-7 




Q4 ~i~ Q4 


-2 - i 
^ 2 


QlQl 






QlQl 


QIQ7 + Q7QI 


-8-1 



Table 1: Low- lying quasi-primary operators. T is the stress-energy tensor. 



3.2 Computation of the Scaling Functions 

We are now in a position to compute the singular part of the scaling functions of ther- 
modynamic observables, e.g. the free energy, making use of the lattice Hamiltonian we 
have constructed. This is achieved by starting with the partition function Z{he,P) of the 
lattice Hamiltonian Hiat- Expanding this expression, it is possible to write down a formal 
series expansion in the variables hi and ^ ^ 0. Hence we are able to obtain the scaling 
functions expressions for the non-scaling corrections to the observables of Section |1.1| . 
To obtain these results, some remarks are in order: 

• In the (numerical) transfer matrix analysis we are interested in the limit ^ ^ 
(the thermal perturbation is smaller than the magnetic one), so we can consider the 
following expressions for the VEV of the operators 

, . , , dim O , , . . > 

(O) = Ao h-d^ qo{^,he) (49) 
^We report in appendix ^ the expression of some of the scahng functions. 
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where qo is an analytic function of its arguments. In the general case this ansatz is 
not correct because the presence of resonances induce the appearance of logarithms 
in the expression of VEVs. 

• In the formal expansion of the partition function Z{h£,f3), we find that there also 
appear products of the conformal fields contained in Hiat- The correct way to deal 
with these is to use the fusion rules of the conformal theory 

WM-W] (50) 
[a][a]~[l] + [6] 

where the notation [. . .] means that we are referring to the whole conformal family. 



3.2.1 A Peculiar Case: The Free Energy 

In this Section we show how to determine the scaling functions of the model in the region 
of interest ^ 0). We have computed the scaling functions by means of two different 
methods: the renormalization group, and CFT approach. 

If we consider the free energy as an example, in the renormalization group approach, we 
can write it as the sum of three contributions 

fit, he) = fb{t, he) + fsingit, he) + fiogit, he) (51) 

The bulk term takes into account analytic contributions in the variables t and he due to 
non critical behavior 

Mt, he) = C + /o,2 h'e + /o% ht + (/^o + /i,2 + /m ht) t+ 

+ iflo + /2,2 h'e + fl, ht) + ifl, + fl, h] + fl, ht) + 0{hl t'). 

Only even powers of he appear because the free energy is even under Z2 transformations. 
The non-analytic contribution is given by the master equation of the RG 

U.,it, he) = gl^^Y (j^, {g^ ^fl/^} j (53) 

where the scaling variables gh{t, he), gt{t, he), gu{t, he) are defined in the usual way |TD[ 

gt{t, he) = t + bthj + ctt^ + dtt^ + etthj + ftt^ + gthj + htt^hj + ... 

gh{t, he) = he{l + Cht + dht" + Chh] + f^t^ + guth]) + ... (54) 

guit, he) =u + aut + buhj + cj^ + dj'^ + euth] + fuh'l + ij^h] + . . . 

and some of the coefficients are known either exactly or numerically (see appendix |^). 
It has been known for a long time p3[ that, in order to take into account the logarithmic 
divergence of the specific heat of the Onsager solution the term fiog ( pT| ) must have 
the following form 

= gl{t, he) log {\g,{t, he)\-') Y (j^, {g^ g^^^^}^ . (55) 
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The function Y{-,-) can be considered as a constant in our analysis, whose exact value is 

n,-)^ A, = 

The previous expression needs some clarification in order to extend it around the magnetic 
axis. Following the discussion of p3[, we point out that if we want f{t, hi) to be an analytic 
function of t at fixed hi, a logarithmic term should appear in the expansion of Y{-, ■). In 
order to make the log \gt\ term disappear, leaving us with the correct log \gh\, we have to 
extract a term like 

x'^\gh\^ log\x\~^Af (56) 
from the expansion of Y{-, ■) when x — > 0, where x is the adimensional variable -tja- 
Hence the logarithmic contribution has the following form 

f,^g{t,he) = g^\og\g,\-'^. (57) 



This gives the correct divergence of the specific heat pS 



Cr^^r^2^\Og\h,\+.... (58) 

With regard to the CFT approach one remark is in order: the CFT derivation can 
generate only the singular term (/sj„g) , while the bulk and the logarithmic contributions 
must be added separately. 

The comparison of the two different approaches (CFT and Renormalization Group) 
enables us to understand the relation between the operators and the corrections they give 
rise to (for a careful discussion about this crucial point see 0). A relevant example is 
given by the stress-energy tensor: it is the only responsible for terms like /i^^/^^, t h?"^/^^ , 
^2 /^32/i5 fjj-^jg have fundamental implications in discussing the numerical data. 



4 Analysis of Results 

The analysis of the data obtained by transfer matrix technique is done as follows: 

• To analyze the large amount of data, due to the large number of required terms in 
the expansion of the scaling functions, and in order to obtain sensible results in the 
fits, we fix as many terms as possible, resorting to both exact and high precision 
numerical results already known in the literature; 

• Fixing all known parameters in the scaling function and then fitting the data enables 
us to check their correctness in this particular regime (^ — ^ 0, i.e. the neighborhood 
of the magnetic axis); 

• The high precision of our data enables us to conjecture the presence or the absence 
of contributions due to well identified sources, i.e. terms like {At-\-Bt^) h^"^^^^, which 
are entirely due to the stress-energy tensor (the irrelevant fields TT, T^ and T). 

This Section is devoted to developing such an analysis and to discuss our results. 

•^In principle it is possible to include higher logarithmic powers in the free energy scaling function, 
However the leading order of the expansion of such terms are zero with a precision of 10^^ in our fits. 
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t 


he 


L 


0.0130481 


0.01 


9 


0.0119726 


0.02 


10 


0.0108971 


0.03 


11 


0.0087367 


0.04 


12 


U.UU r OOio 






0.0065668 






0.0043875 


0.17 


18 


0.0021985 


0.18 


19 


0.0010992 


0.19 


20 


0.0 


0.20 


21 



Table 2: Parameters value. 



4.1 Outline of the Numerical Computations 

We perform our numerical computations extracting the four larger eigenvalues of the 
transfer matrix at given t, he, and L, the size of the transverse direction of the lattice. 



The diagonalization method utilizes an iterative algorithm, due to that evaluates the 
highest-lying eigenvalues of the transfer matrix with arbitrary high precision. The choice 
of an iterative algorithm lies in the reason that the typical size of the reduced transfer 
matrix is of order 10^ and the exact diagonalization, even if partial, is an unobtainable 
computationally. 

Our programs we run on 10 PCs equipped with Pentium III processors and 256 Mb of 
RAM for a total CPU time of about 2 months. 

We performed our numerical computations for each available choice of the three parame- 
ters t, he and L that we collected in table ^ for a total of 2600 different runs. 

4.2 Infinite Volume Extrapolations 

After obtaining the values of the thermodynamical quantities at fixed L, we would like 
to extrapolate them, taking the thermodynamic limit. In order to do this, we follow the 
method of JlOl, which we outline briefly for sake of completeness. 



The basic idea is that the behavior of any thermodynamic quantity in a massive QFT will 
show an exponential decay as a function of L. The task we would like to achieve is to find 
the asymptotic value of the observables; in order to reach this result, we must be able to 
subtract all the possible exponential behaviors from our data. To achieve we iteratively 
subtract the exponential behavior from our data, removing both leading and subleading 
corrections. We define 



b\Li - 2) = c exp(-x(Li - 2)) + 

b'{Li - 1) = c exp(-x(Li - 1)) + b'+\Li) (59) 
b\L^) = c exp(-x(Li)) + b'+\L^) 
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where b^{Li) is the transfer matrix quantity. A step of iteration is defined by solving 
the previous system with respect to c and x. The iteration chain stops when 

the predicted value becomes numerically unstable. The taken result is the last stable 
prediction and its error is evaluated from the variation with respect to previous step. 

4.3 The Fitting Procedure 

In order to analyze our data we fit them with the scaling function previously obtained for 

a generic operator O. 

The fit is performed in three steps: 

1. We start fitting the data at fixed he with a polynomial in the reduced temperature 

0{t)q\h, - O(0)U, = Afihe) t + A^{h,) e + A'iih,) + (/i,) t^ + ... (60) 

where we found convenient to subtract the value of the observable at the critical 
temperature; with this subtraction we are a position to discard all the terms in the 
expansion that depend on only. The results of the fit are collected in a table 
displaying the values of Aiihi) and ^^(/i^) as a functions of h^. 

2. By means of the expansion for the scaling function, we can fit the functions Aiihi) 
and A2{he) against the magnetic field h^- In this way we were able to obtain the 
amplitudes of the corrections to the observable O at first and second order in the 
reduced temperature t. 

3. We use all predicted values of the parameters of the fit to calculate the with a 
unique fitting function for all the data (the coefficients A^{he) and A^{he) of step 
1 are only useful to the determination of the of the entire set of data). 

The list of requirements a fit must fulfill to be accepted is the following: 

• The reduced of the fit had to be of order unity, i.e. we required it to have a 
confidence level larger than 30%. 

• In order to be included in the fitting function, the subleading terms must have an 
amplitude larger than the corresponding error. 

• The number of degrees of freedom must be larger than 5. 

The previous requirements are very hard to fulfill at the same time; this way we also take 
into account the systematic errors due to extrapolations. 

4.4 Numerical Results 
4.4.1 Determination of Ch 

The least precise of the constants collected in the appendix ^ is Ch- The high precision 
of our data enables us to make an attempt to improve this estimate. Following the same 
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procedure of , and fixing all the known terms in the scaling function we are able to give 
the following result 

e^^ = -0.007298(3) (61) 
where the data we used as input are collected in appendix ^ 

4.4.2 Free Energy 

The analysis of the free energy can be performed in four steps. 
Firstly, we fit the the term proportional to t with the following function 



_8_ 16 

- 0.623226 - 0.511645 hi"- + 0.329993 hV' - 0.0427280 



- 0.0107535 h] log + h}^ Kq + 0.0184589 /if + h^' /Ci+ (62) 

46 48 52 

+ 0.0001605 h^" + hf /C2 + hf /C3 + . . . . 

As we explain in the appendix ^ we set y'^^'^^(0,0) = according to the results 
of |2^ and moreover we can found y*^^'^-'(0, 0) compatible with zero within the 



error in all our fits. 

The fit shows a non-zero correction due to the stress-energy tensor, in fact we are 
able to estimate the value of /Co 

- 0.05828 < /Co < -0.05820. (63) 

where the expression of /Co in terms of the scaling fields is given by 

/Co = (a„ + 0.664773m) ^("'^^(O, 0). (64) 

This is quite an interesting result, in fact it is known that in both the integrable 
perturbation of the critical Ising model, the amplitude of this term is compatible 
with zero. 

Secondly, we perform the fit of term with the following function 

_8_ 16 

/C4 - 0.131877 \og\hi\ - 0.244467 hf + 0.299064 hf + 

22 32 V^^J 

+ /C5 hf - 0.0318458 h] + 0.0100527 h] \og\he] + hf + . . . . 

As previously stated we put F(^'^)(0, 0) = F(^'^)(0, 0) = 0, and we found that leading 
correction is 

0.03655 < /C4 < 0.03657. (66) 
The expression of /C4 in terms of the scaling fields is: 

/C4 = 0.481290866 + f^^ (67) 
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in this way we are able to give an estimate for the bulk term 



_ 0.44474 < f^Q < -0.44472. (68) 

Furthermore, we can clearly see that the term due to the stress-energy tensor )Cq 
(which comes with the power h^^^^^) is different from zero, as for the previous case. 
As discussed in it is difficult to have reliable estimates for the amplitudes of the 
subleading termsQ due to large systematic deviations induced by the uncertainties in 
the leading corrections, nevertheless we can still assert that they are different from 
zero. 

Thirdly, we check that within the precision of our computation, the function Y{-,-) is a 
constant. In fact, if we expand it in Taylor series, we find that the first contribution 
to the scaling function comes in the term and has the following form 

F(^'°)(0,0) \og\h,\ hf''\ (70) 

If we perform the fit with this new contribution, we find that its amplitude is 
compatible with zero. Hence, for our purposes it is safe to consider Y{-,-) as a 
constant. 



Fourthly, we compute the x with the global scaling function (Eq. ^O]) both on t and /i^ 
utilizing the value of the constants predicted in the previous steps, in order to verify 
the correctness of the fitting procedure. 

It is important to mention that in the fitting functions we use the known values of the 
constants of Appendix 0. It is a non-trivial test on the validity of both our results, and 
the well known values we used. 

4.4.3 Magnetization 

Firstly, we write down the correction proportional to the t term 

2728775 J- 
- + 0.3519930 hj" - 0.0747026 + 0.0215069 h^log \he\ + 



h 
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17 23 5 ai 11 f71") 

+ hY^Ho + 0.0467625 + hJUi - 0.0004923 + V ^^2+ 

37 41 

+ hfn^ + hYni + ... 

where all the known quantities are taken into account. In particular we also im- 
pose the constraints we found in the analysis of the free energy, i.e. ^'•^'^''(O, 0) = 



"^As discussed in jg], notwithstanding systematic errors, one can give also a rough estimate of the 
subleading amplitudes. Following the same route we are able to find 

- 0.72f < /Cg < -0.717. (69) 
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Y^^'^\0, 0) = 0, and we find perfect agreement also in tliis case. 

Tlie leading correction is due to the stress-energy tensor, as we expected from the 

analysis of the free energy, and its amplitude is given by 

- 0.124 < 7^o < -0.122 (72) 
which is compatible with Tig = f| within our errors. 

Secondly, the contribution due to t"^ is given by 

0.1318769 0.1303825 i 22 ^ 
^ + 0.3190016 h^' + — /if - 0.0536389 he+ 

+ 0.0201054 hi log \ht\ + hf He + . . . 

and is obtained following the same strategy as before. This order also shows a non- 
zero contribution due to stress-energy tensor, i.e. Hq. 
The amplitude of leading order correction is given by 

- 0.723 < 7^5 < -0.718. (74) 

We can check that it is in reasonable agreement with the subleading amplitude of 
the free energy = /C5 (see (pD|)). 

4.4.4 Internal Energy 

Firstly, the scaling function for the t term is 

_8_ 16 22 

Wo + 0.2992532 log \he\ + 0.5547414 + hf^Wi + /if W2 - 0.0918396 /i^+ 

32 38 

+ /if W3 + /if W4 - 0.0228115 h'j log \he\ + ... 

the predicted value for the leading contribution is 

- 0.1659 < Wo < -0.1657 (76) 

which expressed in term of scaling function is Wq = (—2.1842763 + 4.5383706 /|o)) 
the value of /20 is consistent with the previous estimate (^Sj). As before the fit is 
compatible with the constraint ¥^'^•'^^0,0) = F(^'^)(0,0) = 0. 

Secondly, the scaling function for the t"^ term is 

W'^ A M 16 

^ + We - 0.2797532 log \he\ + /if W7 + /if We + /if Ws+ 

/if (77) 

22 

-0.0727482 /if W5 + . . . 
the predicted value for the leading contribution is 

- 0.414 < W5 < -0.408 (78) 

where W5 = 0.5672963 Y^^'°\0,0), and again the fit is compatible with the con- 
straint r(2,i)(o,o) = r(i'i)(o,o) = o. 
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4.4.5 Susceptibility 



The form of the scahng function is in total agreement with the computed values of the 
this observable. However no new predictions are available. 



4.4.6 Mass spectrum 

In order to check the results of the FFPT proposed by we performed the fit on the 
square of the first three masses of the theory, and computed the following ratio 

5ml m(t)f-m(0)f ^ ' 

where m{t)f are the perturbed masses, and m{0)f the unperturbed ones. 
The results of our analysis are reported below: 

= 1-393(17) 



''3 



Smf 



3.16(30). 



These results are in perfect agreement with the theoretical prediction. In fact, since the 
deviation from the integrable model is small, it is correct to write 6m'^ = 2m{0)i6mi, and 
finally we have 

P = 0.86(1) 

— ^ = 1.58 15). 

omi 

4.4.7 Magnetization Overlap 

The magnetic overlap P can also be analyzed as before. We are able to find the leading 
order corrections of terms proportional to t 

Tin , 2 „ 

— ^ + h^Ui + /^T57^2 + /i^7^3 + /lw7^4 + h^^n^ + /^57^6 + hTBn^ + h^Us + ... (82) 



and 



TZq 7?-in ^^^^ , 2„ , 8 „ , i4„ , > 

-^ + ^ + ^ + hT^TZu + h->ni3 + h^TZu + • • • • (83) 

hl5 hl5 /il5 

The numerical values are 

0.628 < 7^o < 0.631 , , 

84 

0.661 < 7^9 < 0.664 ^ ^ 

It would be interesting to calculate the same corrections on theoretical grounds (at least 
for the t correction) in order to make a comparison. 
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5 Discussion and Conclusions 



In this paper we studied the effect of a mixed relevant perturbation on the Ising model 
using the Transfer Matrix technique. We consider the neighborhood of the magnetic axis 
in the limit ^ -C 1, where ^ is the adimensional parameter t/h^^^^. 
We have concentrated our efforts on the following 

• We calculated the scaling functions (appendix |BD applying a CFT approach to the 
problem. 

• We used all known predictions about the behavior of scaling functions (see appendix 

and verified that they all agree with our data. 

• Being able to identify the contribution of secondary fields to the scaling function, we 
predicted that there is a non zero contribution due to the stress-energy tensor, and 
we evaluate it. In our opinion this a quite interesting result, because it is known 
that the contribution of these particular secondary fields is zero if we study the 
model with only one relevant perturbation. 

• We obtained estimates of several amplitudes never predicted by any other analytical 
method before. 

• We calculated the correction to the Zamolodchikov mass spectrum of the Ising model 



with a magnetic field, and we found perfect agreement with Delfino et al. 14 



There are two possible developments of this work: We can use our data to improve the 
knowledge of the equation of state of the Ising model, in order to map all the possible 
regime of perturbations. We can extend the analysis of ^ to study the effect of mixed 
perturbation on the finite temperature results. 
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A Known Numbers 



It is known that it is possible to write the free energy of the model in terms of nonlinear 
scaling fields . 

The scaling fields are analytic functions of t and he respecting the Z2 parity of h^. 
Their Taylor expansions are expected to be 

gt{t, he) =t + hth] + ctt^ + dtt^ + ctth] + /tt^ + c/t/i^ + htt^h] + ... 

Qhit, he) = he{l + Cht + dht" + e^hl + f^t^ + ghth1)+ (85) 

gu{t, hi) =u + aj + huh] + Cyf + dut^ + ej,h] + + l-ufh] + ... 



Here we report all the analytically known coefficients |26, 2S,29, 30 



_ /3c , 23/52 191/?f 



72' " 16 ' 48x/2' 

c - d-^^' f- (86) 

^ ^ 



6^/3,72, 



'16/32 



where 



= 0.0403255003 . . . (87) 



From the analysis of the model at the critical temperature and h^ ^ 0, we obtain a new 
estimate for eh 

eh = -0.007298(3). (88) 

For the free energy, we use also the coefficient reported in the following 

F(0,0) = 0.99279949... 
y(i'0)(0,0) = A^f^Jbt = -0.511645336... 

A\ 2 = , ^ = 0.0208602 . . . (oq\ 
^'^ F(0,0)8/3c y^^! 

/,= -G + ^log2 

TT 2 

A) ^^ = -0.0524442 . . . 
where G is the Catalan constant and the A's are defined in f^. 

A high precision study of the thermal perturbation have been performed by Orrick et al. 
in p6|], from this work we are able to extract a set of parameters for our scaling functions. 
In particular we are able to observe that the known coefficients for both thermal and 
magnetic perturbations are exactly consistent, and we can extend the knowledge about 
the thermal coefficients to gain predictions about the unknown magnetic terms. 
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To achieve these results we compare the scahng function of the susceptibihty along the 
thermal axis with the high precision estimates of pGf and pSf, in this way we obtained 



y(2.o)(o,0) = 0.9625817322. . . 

(90) 

/2% = 0.000528775 ... . 



/o,2 = ^(0, 0) = -0.05206662255 . . . 

f^^ = -0.00348278... 



These results are obtained for the thermal perturbation in the regime t ^ 0, hi = 0, we 
observe that these results are valid also in our regime of interest, because the analytic 
continuation of Section ^.2.1| do not affect similar kinds of terms, in total agreement with 
our fits. 

Furthermore we obtain also the relations: 



F^^'^) (0, 0)(a„ + 0.701128m) = 
y(2.i) (o,0)m = 



(91) 



In order to fulfill the above requirements, we set F (^'^^ (0, 0) = because the other choice 
u = 0, ttu = was not consistent with our fits. 

Finally, from Onsager's exact solution |jl]] we shall extract 

/{-Q = -0.623226 .... (92) 



B Scaling Functions 

Here we report some of the scaling functions we obtained. We remark that the functions 
have been evaluated to higher order of the expansion reported here (0(/i|)0(t^)). 
Free Energy: 

fit)\h, - /(0)U, = Aiih) t + Aiih) t' + ... (93) 



_8_ Ifi 16 8 

+ (^-^ Af log \h,\ h + fl^ + h Y'^'''\0, 0)^ h^^ + ... 



(94) 
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1 22 / 8 \ 8 

+ -bt y (0, 0) /if + f au y (0, 0) + - II Cft y (I'l) (0, 0)+uct Y^^''\0, 0) j /i| + 

(16 16 8 16 

-Y^AfbtCh- Y^Af log \hi\ btct- —AfCh- — Af log \he\ et + /2_2 + 



+ (0,0) + 6^^(2,0) (0,0)) /i; 



(95) 



Internal Energy: 

E{t)\h^ - E{0)\n, = Afik) t + A^{h) e + (96) 
(M = (-^ ^/ log M To + 2 flo Tc + Tey(^'°)(0, 0)) + 

+ (^^c,r,r(^'0)(o,o) + 2Qr,r(i'°)(o,o)^ hp+ 

(1 ry '^0 \ 16 

— c/,2Ter(0,0) + -4TeF(0,0)+«TeF(^'i)(0,0)J hf+ (97) 

+ bt T, (0, 0) hf + i2 a, F^^'^) (0, 0) + 

+ ^uchT, y (^'1) (0, 0) + 2uctT, y (0, 0)^4 + ... 
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A^ih) = -\ Y^'''^ (0, 0)+(-^Ajc,n-^Af log \he\ CtT, + 3 f^, T,+ 
2 /if V ^ ^ 

+ 3c,T,F(2'°)(0,0)) + (^-|c,^T,F(i'°)(0,0) + ^c,QT,r(^'°)(0,0)+ 

+ htnY^'''\o, 0) /if + (-^ c,3T,y(o, 0) + ^ c, 47;y(o, o)+ 

+ -fhT, y (0, 0) + - a„ Te y (2,1) (0, 0)+ uCkT, y (^-i) (0, 0)+ 
o z o 

+ 3«QT,y(2'i)(0,0)j /if + (^-^6iC;,T,y(3'0)(0,0) + 3 6iQT,y(3.o)(0,0) + 

- le,T,y(3'°)(0,0) + ^e,T,y(-^'°)(0,0)^ /,i + a„ y(^'^)(0, 0) + 

+ ^ c^' y ^'-'^ (0, 0) + 3 a„ ct Te y ^'''^ (0, 0) + ^ ii c* 7; y ^'''^ (0, 0)+ 

24 

+ 3 c„ y (1'^) (0, 0) + — M 4 Te y (0, 0) + 3udtT, y (1'^) (0, 0) + 

o 

+ ^ii'7;y(3'2)(o,o)^ /i| + ... 

(98) 

C Example of numerical data 

Here we report, for sake of completeness, an example of the data we used to fit. 
We refer to value oi p = 0.4373. 
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fh Ejf.lu) 

0.01 0.035901(3) 

0.02 0.0569877(2) 

0.03 0.07237192(2) 

0.04 0.084838424(2) 

0.05 0.0954482541(1) 

0.06 0.1047420407(1) 

0.07 0.11303997732(1) 

0.08 0.12055040575(1) 

0.09 0.12741797596(1) 

0.10 0.13374801551(1) 

0.11 0.13962007086(2) 

0.12 0.145095979793(2) 

0.13 0.150224955365(1) 

0.14 0.155046932292(1) 

0.15 0.159594851224(1) 

0.16 0.16389626622(1) 

0.17 0.167974505628(1) 

0.18 0.171849529388(1) 

0.19 0.17553857454(1) 

0.20 0.179056649653(1) 



0.931944971(5) 

0.939788377(1) 

0.94794151643(2 

0.95629353642(1 

0.96479217393(1 

0.97340669611(1 

0.98211679357(1 

0.99090802119(1 

0.999769568339(1 

1.008693035843(1 

1.017671707919(1 

1.026700091153(1 

1.035773608172(1 

1.044888386026(1 

1.054041105248(1 

1.06322888925(1) 

1.072449221416(1 

1.081699881718(1 

1.090978897417(1 

1.100284504134(1 



0.76120(4) 

0.802572(1) 

0.82640939(1) 

0.843168904(1) 

0.8560620814(1) 

0.8665107103(1) 

0.8752714873(1) 

0.8827957916(1) 

0.88937481546(1) 

0.89520749914(1) 

0.90043585646(1) 

0.905164845833(3) 

0.909474255678(2) 

0.913426172207(2) 

0.917069865355(1) 

0.920445095734(1) 

0.923584418311(1) 

0.926514827477(1) 

0.92925895744(1) 

0.931835974906(1) 



1.45013(2) 
1.150734(1) 

0.9510868(1) 

0.808372(1) 

0.701236(1) 

0.6178284(1) 

0.5510433(1) 

0.4963581(1) 

0.4507573(1) 

0.4121515(1) 

0.3790476(1) 

0.35034994(1) 

0.32523603(1) 

0.303076426(1) 

0.283381176(1) 

0.2657631501(1) 

0.2499123945(1) 



hi 


l/mi(i, he) 


l/m2{t, hi) 


l/m3(i, he) 


Frit, he) ' 


0.03 


1.66609(4) 








0.04 


1.42446(2) 


0.86(2) 




0.413(1) 


0.05 


1.262217(1) 


0.775(2) 




0.403172(4) 


0.06 


1.1438966(2) 


0.704762(5) 




0.39427(1) 


0.07 


1.0528411(3) 


0.64979(3) 


0.5334(3) 


0.3859807(2) 


0.08 


0.980049866(1) 


0.60587(2) 


0.4973(3) 


0.37813145(1) 


0.09 


0.920186381(1) 


0.569787(5) 


0.4676(3) 


0.3706210923(3) 


0.10 


0.8698623782(3) 


0.539467(1) 


0.4417(3) 


0.3633856671(3) 


0.11 


0.82681092694(1) 


0.5135465(5) 


0.4203(2) 


0.356382554(1) 


0.12 


0.789451203443(1) 


0.4910712(1) 


0.4019(1) 


0.3495820201(7) 


0.13 


0.756643166957(1) 


0.47135216(1) 


0.38589(7) 


0.34296250455(5) 


0.14 


0.727541458508(3) 


0.453877992(1) 


0.37178(2) 


0.33650783473(4) 


0.15 


0.701504369616(1) 


0.4382602553(3) 


0.35921(1) 


0.3302055233(2) 


0.16 


0.678034937921(1) 


0.424197853(2) 


0.347936(2) 


0.324045673287(8) 


0.17 


0.656741586699(1) 


0.41145329066(3) 


0.3377387(2) 


0.318020257612(3) 


0.18 


0.637311086114(2) 


0.39983637405(1) 


0.3284631(3) 


0.31212262996(1) 


0.19 


0.619489526443(1) 


0.389192731029(1) 


0.31998(1) 


0.306347184601(2) 


0.20 


0.603068643116(1) 


0.379395558407(1) 


0.312184(1) 


0.300689114577(1) 
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